How Well a Chaotic Quantum System 
Can Retain Memory of Its Initial State? 
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In classical mechanics the local exponential instability effaces the memory of initial conditions and 
leads to practical irreversibility. In striking contrast, quantum mechanics appears to exhibit strong 
memory of the initial state. We relate the latter fact to the low (at most linear) rate with which the 
system's Wigner function gets during evolution more and more complicated structure and establish 
existence of a critical strength of external influence below which such a memory still survives. 
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Strong numerical evidence has been obtained [l[ that 
the quantum evolution of classically chaotic systems is 
quite stable, in sharp contrast to the extreme sensitiv- 
ity of classical dynamics to initial conditions and weak 
external perturbations. Being the very essence of classi- 
cal dynamical chaos, this sensitivity results in rapid loss 
of memory and practical irreversibility of classical mo- 
tion. Qualitatively, this crucial difference is explained 
by a much simpler structure of quantum states as com- 
pared to the extraordinary complexity of random and 
unpredictable classical trajectories 0, more rigor- 

ous terms, chaotic classical systems are characterized by 
positive algorithmic complexity described by Lyapunov 
exponent. Unfortunately, being formulated in such a way 
the concept of complexity cannot be transferred, sic et 
simpliciter, to quantum mechanics, where the very no- 
tion of trajectory is irrelevant and there is no quantum 
analogue to the Lyapunov exponent. Therefore, at first 
glance, there exists no quantitative measure of compar- 
ative complexity of classical and quantum states of mo- 
tion 

However, individual classical trajectories are, in 
essence, of minor interest if the motion is chaotic. They 
all are alike in this case and rather behavior of manifolds 
of them carries really valuable information. Therefore the 
methods of the phase space and the Liouvillc form of the 
classical mechanics become the most adequate. It is very 
important that, opposite to the classical trajectories, the 
classical phase space distribution and the Liouville equa- 
tion have direct quantum analogs. Hence, a comparison 
between classical and quantum dynamics can be made 
by studying the evolutions in time of the classical and 
quantum phase space distributions expressed in similar 
canonical variables and both ruled by linear equations. 

The paramount property of the classical dynamical 
chaos is the exponentially fast structuring of the system's 
phase space on finer and finer scales. On the contrary, 
degree of structuring of the corresponding quantum " dis- 
tribution" is restricted by the quantization of the phase 
space. This makes Wigner function relatively "simple" as 



compared to its classical counterpart. The great advan- 
tage of the phase space approach is also that one operates 
with distributions which can be presented in both clas- 
sical as well as quantum cases in identical action - angle 
(/, 9) variables. 

In practice, to explore stability or reversibility of mo- 
tion one computes fidelity 0, 0| which characterizes the 
weighted-mean distance between two distributions evolv- 
ing, for example, under two slightly different Hamilto- 
nians. In this paper, we exploit a somewhat different 
aspect of the Peres fidelity. Namely, we use it to estab- 
lish a quantitative relation between the degree of loss of 
memory on the one hand and complexity of the corre- 
sponding phase space distribution on the other hand. In 
particular, we show how the number M{t) of 9 -Fourier 
harmonics (see below Eqs. ([71[8])) can be used as a suit- 
able measure of complexity of a distribution at a given 
time t. In the case of classical chaotic dynamics, this 
number grows exponentially in time with the rate re- 
lated to the rate of local exponential instability Q. At 
the same time, this number increases only power-like if 
the motion is regular. Thus the rate of growth of the 
number of harmonics is, similar to the Lyapunov expo- 
nent, a characteristic measure of classical complexity. In 
the framework of the phase space approach, the number 
of harmonics of the Wigner function appears, contrary 
to the Lyapunov exponent, to be a relevant measure of 
complexity also in quantum case. In what follows, we ex- 
amine the time behavior of this quantity and its relation 
to fidelity and reversibility properties. A detailed deriva- 
tion of the results reported in this paper can be found in 
Rcf. 0. 

Let H = H{a\a;t) = H^°\n = a^a) + H^^\a\a-t) 
be the Hamiltonian of a generic nonlinear system with a 
bounded below discrete energy spectrum En'' > 0, which 
is driven by a time-dependent force of such a kind that 
the classical motion exhibits a transition from integrable 
to chaotic behavior when the strength of the driving 
force is increased. Here a\a are the bosonic, creation- 
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annihilation operators: [a, a)] = 1 . We use the method 
of c-number a-phase space borrowed from the quantum 
optics (see for example 



M nil) 



It is basically built 

upon the basis of the coherent states \a) = D (^-^^ |0) 

obtained from the ground state |0) of the unperturbed 
Hamiltonian with the help of the unitary displacement 
operator D (A) = exp(Ad^ — X*d). Here a is a complex 
phase space variable independent of the effective Planck's 
constant h. 

The Wigner function W in the a-phasc plane is defined 
by the following Fourier transformation 



W{a\a;t)^ — 



1 



d2^e^''*^~''7l)Tr 



(1) 

where p is the density operator and the integration 
runs over the complex r^-plane. The Wigner function 
is normalized to unity, J (PaW{a* , a;t) — 1 and is real 
though, in general, not positive definite. It satisfies the 
evolution equation 

i^^W{a\a;t)=CgW{a\a;t), (2) 

with the Hermitian "quantum Liouville operator" Cq. 
This equation reduces in the case ?i = to the clas- 
sical Liouville equation with respect to the canonical 
pair a, ia* with the classical Hamiltonian function be- 
ing given by the diagonal matrix elements Hc{a* , a; t) = 
{a\H^-'^\a' , a)\a) of the normal form iJ^^^ of the quan- 
tum Hamiltonian operator. In other words, this function 
is obtained from the quantum Hamiltonian by substitut- 
ing d — > a/v^, ^ a* l\fh. 

We define the harmonic's amplitudes WmiJ) as the 
Fourier components of the Wigner function with respect 
to the angle variable d introduced by the canonical trans- 
formation a. = \fl e~^^ . The normalization condition re- 
duces then to dl Wo(/; ~\ whereas the amplitudes 
of other harmonics are expressed in terms of the matrix 
elements (n -I- •m\p\n] of the density operator along the 
TOth collateral diagonal as 



[U/ny L"^' (4//n) {n + m\p{t)\n), m > 0, 

with L™ Laguerre polynomials and W—m = W^. The 
inverted relation reads 



(n + m|pW|n) = (-ir2yg=;x 



(4) 



Aiming to connect the reversibility of the motion with 
the complexity of the Wigner function, we follow the ap- 
proach developed in Ref. We consider first the for- 
ward evolution p{t) = U{t)p{0)W{t) of a simple initial 
(generally mixed) state p{0) up to some time t ~ T. 



An instantaneous Hermitian perturbation S,V S(t — T) 
with the intensity ^ is then applied which transforms the 
state /3(r) into yo(T,^) = P{C)p{T)P^^) . The resulting 



transformation P(^) 



is unitary. For example 



this transformation is equivalent to the global rotation 
W{I, 0; T) W{1, 9 + C;T) by the angle ^ in the phase 
plane if the operator V ~ n. In particular, we use be- 
low an infinitesimal perturbation of such a kind to reveal 
complexity of the Wigner function at the instant T. 

The new state p{T, ^) serves as the initial condition for 
the backward evolution U{—T) = W{T) during the same 
time T, after which the reversed state 

P{0\T,0 - (fHT)KT,OU{T) ^ P(e,T)p(0)P^(e,T), 

(5) 

is finally obtained. Here P{£,,T) = e~*«^(^\ with 
V{T) = U^T)VU{T) being the Heisenberg evolution of 
the perturbation during the time T. At last, we con- 
sider the distance between the reversed /5(0|r, ^) and 
the initial p{Q) states, as measured by the Peres fidelity 
Frev{^;T) = Tr[;5(0|r,e)p(0)]/Tr[p2(0)] i which can be 
also expressed with the help of Eq. ^ in terms of the 
Wigner function as 0, 



^rev (C? ^) 



/ d^aW{a*,a;0)W{a\a;0\T, £_) 



J d^aW^ {a*,a;Q) 
J d^a W {a* ,a;T)W {a* , a; T, C) 



(6) 

The fidelity is bounded in the interval [0, 1] and is the 
closer to unity the more similar are the initial and re- 
versed states. At the second line, which is a consequence 
of the unitary time evolution, the fidelity F{^;T) mea- 
sures the complexity of the Wigner function at the mo- 
ment t — T (see below). Both the functions Frev{^',T) 
and F{^;T) are numerically identical. The relation ([6]) 
plays the key role in our analyses. It allows us not only 
to relate the degree of reversibility to the complexity of 
the state at the reversal time T but also to establish a 
strong restriction on the upgrowth of the number M{t) 
of harmonics of the Wigner function. The crucial point 
is that while in classical mechanics the number of Fourier 
components has no direct physical meaning, in quantum 
mechanics the number of the components of the Wigner 
function at any given time is related to the degree of exci- 
tation of the system (see for example eq. (jl2p below) and 
therefore unrestricted exponential growth of this num- 
ber is not allowed [1, [H, [3l • It is important that the 
representation ([6]) is valid not only in the quantum but 
also in the classical case when the Wigner function re- 
duces to the classical phase-space distribution function 
Wc{a*,a;t). 

First of all, let us establish a connection between fi- 
delity and complexity of the Wigner function at an ar- 
bitrary moment t. To do this we make the instant rota- 
tion e~*^" at this moment and utilize then the expression 
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given in the second line in Eq. ([6]). This yields Q 
Fie, = 1-2 sin^ (em/2) = 



or): 



where 



Ei=orn^W„M and 



Wm>o{t) 



{2-5rao)j;^dI\Wm{I;t)\^' 



E 



+ 00 



(7) 



(8) 



We define the number of developed by the time t har- 
monics of the Wigner function of a classically chaotic 
quantum system as A4 (t) = ^ (m^)^ what is in line with 
the refs. [1, The set of positive definite quantities 
>V„i is normalized to unity, Em'^o ^™ — ^^"^ '^^'^ 
therefore given a probabilistic interpretation. Notice that 
our definition of the number M. (€) is applicable to clas- 
sical dynamics as well, provided that the harmonics of 
the classical distribution function Wc instead of those of 
the Wigner function W are used in Eq. fH]). Hence, the 
number of harmonics at arbitrary time t is defined via 
Peres fidelity as 



de 



(9) 



As an illustrative example we consider further the 
kicked quartic oscillator defined by the Hamiltonian 



H{a\a) = h 



ujon - 



t2 -2 

2 n 



Vhgit){c 



(10) 



where g{t) = go Es ^(^ ~ ■^)- our units, the time and 
parameters h, luq as well as the strength of the driving 
force are dimensionless. The period of the driving force 
is set to one. The classical dynamics of such an oscillator 
becomes chaotic when the kicking strength go exceeds a 
critical value 50, c ~ 1- The angular phase correlations de- 
cay in this case exponentially and the mean action grows 
diffusively with the diffusion coefficient D ^ g^. 

In the chosen model, the only difference be- 
tween the classical and corresponding quantum 
Liouville operators consists in the substitution 



Wo + 2|a|' 



UJo 



2 da* da 



+ 2\a\ 



{a*d/da* — ad/da) in the unperturbed {go ~ 0) part 
L^^^ . Opposite to the continuous hmction 2|ap, the 

spectrum of the operator K ~ — ^fi^ da* da ~^ ^l*^!^' 
which coincides formally with the Hamiltonian of a 
two-dimensional linear oscillator, is discrete. This 
reflects the quantization of the phase space which stays 
ultima analyst behind the much slower growth of the 
number of Harmonics of the Wigner function than that 
of the corresponding classical phase space distribution. 

A numerical illustration of this statement is given in 
Fig. [T] The initial state is chosen to be a pure ground 
state /5(0) = |0)(0| which corresponds to the isotropic 
Wigner function Ty(a*,Q;;0) = 2e~^l"l with the size 
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FIG. 1: (color online) Root-mean-square {m^)t versus time t 
at go — 1.5. Squares, diamonds and triangles correspond to 
h = 0.01, 0.1 and 1. Empty circles refer to classical dynamics 
and the dashed line fits these data. 



1/2. This size is kept constant throughout all calcula- 
tions whereas the quantum Liouville equation is solved 
for a number of decreasing values of the effective Planck's 
constant h thus approaching the classical dynamics. Let 
us notice that for /i < 1 our initial condition corresponds 
to an incoherent mixture of eigenstates of the operator 
H^^\ It is clearly seen that the exponential increase of 
(TO^)t takes place only up to the Ehrenfest time scale 



I, consistently with the findings reported in 
Then, a much slower power-law increase 



t B cx In h 111 
Refs. d [3 
follows. 

To explain such a behavior we turn to the equiva- 
lent representation of fidelity given in the first line of 
Eq. ([6]). We restrict ourselves to the case of the pure 
ground initial state. More general consideration which 
includes also arbitrary incoherent initial mixtures is pre- 
sented in 0. Using this representation we find in this 
case that Frevi^t) = |/(C,OP where Q 



1 0) = Tr [e-*«"p(t)] 



(11) 



cos(5/2) Jo ^ , tj ■ 

In such a way we relate the behavior of fidelity to evolu- 
tion of the action variable. Comparing now the terms 
in the expansions of the both possible representations ([7|) 
and (jlip of fidelity we arrive at the following significant 
exact relation between the number of harmonics and the 
root-mean-square deviation of the action 

1 



2X2(0, X2it) 



ft2 



{{1% - (/)2) . (12) 



Thorough numerical study Q convince us that af- 
ter proper averaging over strong irregular fluctuations 
(coarse graining) the zero harmonic amplitude decays ex- 
ponentially with the action / at any given time t > tE, 



WoiF,t) 



{I)t + 



(13) 
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FIG. 2: Mean value {I)t/go as a function of time t. Squares 
and triangles correspond to {h, go) = {l,2) and (2,3). The 
straight line shows the classical diffusion law (J)t = got. 
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FIG. 3: {m^)t vs X2(t) ~ ^ + 1 ) , for t = 1, . . . , 100. 

Data correspond to ft = 1 and go = 2. 
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It follows from Eq. ^ that X2{t) ~ % (" 

therefore, M.{t) « The time dependence of the 

mean action (/)< (the deterministic quantum diffusion) 
is shown in Fig. [2] whereas the validity of the stated con- 
nection between the number of the harmonics and the 
excitation of the system is illustrated in Fig[3l Thus the 
number of harmonics grows after the Ehrenfest time not 
faster than linearly. 

It can be readily shown now that for any finite ^ ^ 1 
fidelity equals in the approximation (flS]) 



Frev{i;T)=F{^-t = T) = 



1 



['m?)T 



(14) 



More than that, this formula is valid for any time in- 
cluding the times T < Te 9}. The found result 
shows that a crossover takes place near the critical value 
^c(T) = V2/M{T), from good, F{^;T) « 1, to broken, 
F{^;t) ~ (C/'fc(T))2 ^ reversibility. Our numerical 
simulations (see Fig. H]) nicely confirm the formula (|14[) . 

The dependence of £,c{T) on the reversal time T is 
strikingly different within and outside the Ehrenfest time 
scale. When in the former case it drops exponentially in 
accordance with the classical exponential proliferation of 




FIG. 4: Fidelity F{^;T) versus the scaled variable C/^c(T). 
Data correspond to: (a) h = 1, go = 2; circles: T = 10; trian- 
gles: T = 50. The full curve shows the theoretical prediction 
Eq. (|14|l . The deviations on the tail of the curve are due to 
fluctuations neglected in (|14p . 



the number of harmonics, in the latter case it decreases 
at most linearly. This explains the numerically discov- 
ered [H much weaker sensitivity of quantum dynamics to 
perturbations as compared to classical dynamics. 

To summarize, we have established a quantitative re- 
lation between complexity of the Wigner function and 
degree of reversibility of motion of a classically chaotic 
quantum system. We have analytically proved that the 
number of harmonics M{t) of this function, which can 
serve as a natural measure of the complexity, increases 
after the Ehrenfest time not faster than linearly in strik- 
ing contrast with classical dynamics, where the number of 
harmonics of phase space distribution growths exponen- 
tially. We have shown that if a quantum motion has been 
perturbed at some moment T by an external force with 
intensity ^ and then reversed, its initial state is recovered 
with the accuracy ~ (C/?c)^ as long as the strength is re- 
stricted to the interval < ^ < ^c{T) = V2/M{T). This 
interval decreases with the time T at most linearly be- 
yond the semi-classical domain but shrinks exponentially 
due to the classical exponential instability when the this 
domain is approached. 

The above outlined phase-space approach is quite gen- 
eral and can be readily extended to systems with arbi- 
trary number of degrees of freedom, including qubit sys- 
tems, whose Hamiltonian can be expressed in terms of a 
set of bosonic creation-annihilation operators. Therefore, 
this approach could shed some light of the connection be- 
tween complexity and entanglement, a fundamental issue 
of great relevance for the prospects of quantum informa- 
tion science [1^ . 
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